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Abstract 

In this article we discuss a new Hamiltonian PDE arising from a class 
of equations appearing in the study of magma, partially molten rock, in 
the Earth's interior. Under physically justifiable simplifications, a scalar, 
nonlinear, degenerate, dispersive wave equation may be derived to de- 
scribe the evolution of (j), the fraction of molten rock by volume, in the 
Earth. These equations have two power nonlinearities which specify the 
constitutive realitions for bulk viscosity and permeability in terms of 4>- 
Previously, they have been shown to admit solitary wave solutions. For 
a particular relation between exponents, we observe the equation to be 
Hamiltonian; it can be viewed as a generalization of the Benjamin-Bona- 
Mahoney equation. We prove that the solitary waves are nonlinearly 
stable, by showing that they are constrained local minimizers of an appro- 
priate time-invariant Lyapunov functional. A consequence is an extension 
of the regime of global in time well-posedness for this class of equations 
to (large) data, which include a neighborhood of a solitary wave. Finally, 
we observe that these equations have compactons, solitary traveling waves 
with compact spatial support at each time. 
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1 Introduction 



Consistent, macroscopic models of magma, partially molten rock, in the Earth's 
interior can were developed in [9l[T8l[T9], coupling the flow of the solid rock with 
that of the liquid via conservation of mass and momentum. Under appropriate 
assumptions (small fluid fraction, no large scale shear motions, etc.) the system 
may be reduced to a single scalar equation for the evolution of the fluid fraction, 
the porosity ^, pi [Tl that admit solitary waves. In one spatial dimension, 
the equation is 

5t0 + 5, (0") - 5, {rd. (r'^'M)) = (1.1) 

with the boundary conditions that (j){x,t) — > 1 as a; ±oo. The nonlinearity 
parameter n is specified by a Darcy's Law relationship between the permeability, 
if, of the rock and its porosity of the form K{(j)) oc 0"'. The other nonlinearity 
parameter, m, is related to the bulk viscosity, ^, of the porous rock, with C,{4i) oc 
. It is expected that 2 < n < 5 and < to < 1. A well-posedness theory 
for the initial value problem of (|l.ip is developed in [20]; see also section [5J 

In the article [TT], solitary traveling waves, 4> = 4>c(x — ct), of speed c are 
shown to exist for Ijl.ip for any speed satisfying c > n > 1; to may take any 
real value. In many problems, solitary waves are well-known to be important 
coherent structures, participating in key dynamic processes. In order to play 
this role, solitary waves must be dynamically stable. The most direct approach 
to the nonlinear dynamic stability of solitary waves is via a variational structure 
of the equations. Unfortunately, (|l.ip does not appear to have such a structure 
for the parameter ranges (to, n), arising in the magma problem. However, while 
not of present interest to the problem of magma we observe that when n+m = 0, 
there is a Hamiltonian formulation: 

M = [I -d^ird^ir d,{-r) = j^-^ (1.2) 
m] - 1 + - 1 j (1.3) 

J = J0 = [/-a,(0"9,(0-"-))]~'5„ n + TO = (1.4) 

We will make use of the Hamiltonian structure in these cases to show that 
their solitary waves are orbitally stable^ i.e. for data sufficiently close to a solitary 
wave, the corresponding solution, modulo a time-dependent spatial translation, 
will remain close to the solitary wave in H^(M.). The general method of proof 
is well established and discussed in [21 El HH US |3] for the Kortweg - de Vries 
(KdV), Benjamin - Bona - Mahoney (BBM), and Nonlinear Schrodinger (NLS) 
equations, amongst others. 

We note that in contrast to generalizations of the BBM, KdV, NLS equations 
to arbitrary power nonlinearity, solitary waves of (|1.2p are nonlinear stable for 
arbitrary powers, n. Currently this is established up to a numerical computation 
computation of the slope of the function c i-^ A/'[0c] ; see Proposition [TJ 
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This stability result is also of significance for the global existence theory 
for (|l.ip : at present, no global existence in time result is known for the case 
n + m — 0. We are required to prove, in tandem with the nonlinear stability, 
that solutions for data in a neighborhood of a solitary wave exist for all time. 
Specifically, we note that (|l.ip can potentially become a degenerate dispersive 
equations if (j) tends to zero. As made clear in the well-posedness results [20] 
global existence in time is ensured by uniformly bounding the porosity, (f>, away 
from zero. Solitary waves are examples of solutions, whose porosity is uniformly 
bounded away from zero. The strategy is to show, using spectral and variational 
arguments, that initial data, in a small neighborhood of the solitary wave, 
remain in a small neighborhood, therefore persists in being bounded away from 
zero, ensuring global existence and stability. 

We proceed as follows. In Section [2l we review some of the basic mathemat- 
ical properties of (jl.ip on well-posedness theory and solitary waves, and state 
the our main results: Theorem 12.41 and Corollary [T] In Section [3l we review 
the constants of motion and their relation to the solitary waves. The proofs of 
the main results. Theorem 12.41 and Corollary [TJ on orbital stability and global 
existence of data near a solitary wave solution, are given Section [4] In Section 
[51 we note a relationship between our equations and those that have compacton 
solutions, solitary waves with compact support, [T5], and show that (jl.ip also 
possesses such solutions. 

Finally, we remark that in a forthcoming paper, we prove the asymptotic 
stability of solitary waves in the general case (arbitrary m and n) of (|l.ip . of 
small amplitude, without any restriction on m and n . In fact, the Hamiltonian 
structure which we presently use in the case n + m = for ()l.ip has implications 
for the linear spectral theory and stability analysis in this work, via the Evans 
function (see, for example, [Sj), an analytic function, whose zeros are points 
in the discrete spectrum of and resonances of the linearized spectral problem 
about the solitary wave. 



2 Background and Main Results 

We first state the basic results on well-posedness of (II. ip : see [20], Theorems 
2.12 and 2.13. 

Theorem 2.1 (Local Existence in Time & Continuous Dependence Upon Data) 
Let 00 satisfy 



1 



< 



1 

2^ 



/or i? > 0, 1 > a > 0, and fc > 1. 

(a) There exists Tcxist. > 0, Tcxist. = roxist.(-R, a), and - 1 e C^([0, Toxist.) : 

H'^iM.)), such that cj) is a solution to (jl.ip with t) > a and \\(f>{-, t) — l\\fjk < 
2R for t < Tcxist.. 
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(b) There is a maximal time of existence, Tmax > 0, such that if T^ax < oo, 
then 

lim U{-,t)^l\\H,+ 



(2.1) 



(c) Let 4>^^\ (/i^^' be two solutions of (jl.ip . existing in a common space C^([0, T) 
H^iM.)), T > 0, and satisfying the bounds 



< 2R and 



00) {t) 



1 

< — 

a 



for some < a < 1, i? > 0, fc > 1, j = 1, 2, and all t < T. Then 
there exists a constant, i^cont. — Kcont.{R,Oi,k), such that for any t and 
t',0<t<t'<T, 



< e 



(2.2) 



We will show that the solitary waves of (jl.ip are orbitally stable in the 
following sense. Let us define the distance function, 

Definition 2.2 Let / — 1 and .g — 1 be in H^{M.). Define the sliding metric on 

d{f,g)^M\\f{- + y)-g\\Hi (2.3) 
y 

Definition 2.3 Given (j> - 1 e C\[0,T) : H^{R)), T > 0, we say that cj) is 
orbitally stable, if for all >eO, there exists S > such that if 

||(/.(t = 0)-^(< = 0)11^1 <s 

with^-l e Ci([0,r) : H^R)) 

d((/.(t),V(t)) <e 

for t < T. 

A solitary traveling wave is a solution of the form (j) — (f)c{x — ct), where 
(t>c{x) asymptotes to a constant, say = 1, as a; — s- ±cx). Thus, solitary waves, 
(j)c, of (|l.ip . for the case of n + m = satisfy 

- cdAc + + cd^ {rcdA€d.cl)c)) = (2.4) 

After one integration, and using the boundary condition at oo, 

- c0e + + C(/.;'9,(Ca.(/'c) + c - 1 = (2.5) 

In general, there is no closed form expression for (jjc as a function of a;. As 
previously noted, for c > ti > 1, a solution in excess of the reference state at oo 
exists, (pdx) can be shown to be exponentially decaying as a; — s- ±oo and 
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Figure 1: Two solitary waves for n = 2. Note that as the speed parameter, c, 
increases, the waves become both taher and broader. 



analytic in a strip about the real axis in the complex plane, see |21j for more 

details. Two such waves are pictured in Figure [T] 

It is worth noting, that in [T0l[12], solutions for c = n were found. However, 

(pc{x) = at isolated points and do not fit into our existence theory which relies 

on boundedness away from zero; hence, we do not consider them here. 
(|2.4p also possesses new compacton solutions, discussed in Section [51 
Our Main Theorem and Corollary [1] apply to the Hamiltonian case of 

n + m — 0. 

Theorem 2.4 (Orbital Stability) 

Let 4>c be a solitary wave with c> n and let <f> — 1 G C^([0, T); iJ^(R)) be a 
solution to T > 0. There exists Si, ~ £*{4'c) such that for all e < there 

is a S such that if for some xq G K, 

m;t = 0)-(bci--Xo)\\Hi <S 

then d(0(t),(/)c) < £ for all t G [0,T). 

We make no assumptions about T in this theorem except that T > 0; indeed it 
may be infinite. 

Corollary 1 (Global Existence and Orbital Stability) Given a solitary wave (pc, 
and £ < Si,, there exists 5 > such that if 

1100 - 4>c{- - Xo)\\h^ < S 
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with xo e M, then - 1 G C^([0,oo) : H'^{M.)) and d((/i(t),0c) < £ for all time. 



3 Conservation Laws and Variational Charac- 
terization of Solitary Waves 

3.1 Invariants and Regularity 

In addition to the Hamiltonian Ti, another invariant is the generahzed momen- 
tum 

I + (3.1) 

This was identified in [2D] as a weU defined quantity for (p—lG , formed out 
of a hnear combination of conservation laws discovered in [8j . In Appendix \^ 
we show the relationship between (13. ip and the Lagrangian of 

In order to prove that J\f and Ti. are constant in time for solutions, 
one must first establish their conservation in H^, which is obvious, and then 
approximate an solution in to show that these quantities really are 
invariant. The time Texist. > appearing in Theorem 12.11 is chosen such that 
both Ti and Af are invariant for solutions to (|l.ip . See Sections 4 and 5 of 
pn] for details on the invariance of Af, which can easily be extended to H. 

For (p — 1 e H^, A/" is obviously well defined; H is also well defined; the 
polynomial in the integrand p(a;) — (1— a;""'"^)/(n-|-l)-|-a;— 1, hasp(l) —p'{l) = 
and p"{l) = —n, giving the bound 

\nm<cu-i\\^^ <c 110-111^, 

with C independent of 4>. 



3.2 Variational Characterization of the Solitary Waves 

Let 

For c > n, consider the taylor expansion of £ about a solitary wave (/)c, 



£[^c + u] = £[cP,] + {^f^m,uj + I {^^mu, uj + 0{\\u\\%,) (3.3) 
The variational derivatives are 

= c<j>-r-crdArd:.<j>)-c+i (3.4) 

Lcu = 0[(t)]u = -cd^ {(f>^'^d^u) - [i2n - l)nc02n-2^g^^)2 _^ 2nc(l)^''-^ djc/) + ncf)"-^ ~c]u 

An«2 ^^n„^) _^ _ ^^^n-lg^ (0"9:.0)] U (3.5) 
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[(f)c\ — because of (|2.5p ; a solitary wave of speed c is a critical point of the 
this functional. Alternatively, it can be viewed as critical points of Ti., subject 
to the constraint of Af with Lagrange multiplier c. 

Since the solitary waves are critical points of we would like to be able to 
make an analysis of the form 

\A£\ > i {L,u, u)>C \\u\\%^ + O {\\u\\%) , C> 

to conclude their Lyapunov stability. However, as proved in Proposition [2] , the 
Sturm-Liouville like operator, Lc, is not positive definite; it possesses a negative 
and a zero eigenvalue. Nevertheless, there are natural constraints associated 
with the problem that will ensure positivity; these will be discussed in Section 

For later use, we state a formal result on 



Lemma 3.1 Given a Solitary Wave (j>c, there exist constants D2 and D3 such 
that for all such u £ with < ^• 



£[(j>c + u]- - ^ {LcU, u) 



< D3\\u\\l, (3.6) 



\£[cl), + u]-£[cbc]\ < D2\\u\\l,+D3\\u\\%,=p+{\\u\\^i^.7) 



The polynomial will be used in the proof of the main theorem in Section 

Another property of the invariant J\f evaluated at 0c will imply that solitary 
waves of (|1.1[) are never unstable, see Theorem 12.41 



Proposition 1 (Analytically confirmed forn — 2, suggested numerically n ^ 2) 

J\f[(f>c] > (3.8) 



d 



dc 



for all c > n > 1. 

Proof: Integrating (|2.5p again, and applying the boundary condition at 00, 



2 2 c 



1 - 0g+i 
n + 1 



= f^(0c;c,n) (3.9) 



Using this, and the even symmetry of 

nOC 

N[(t)c] = 2 F{c^c;c,n) + 



ifdx 



(|3.9P can be used to compute 0™'*'^ by solving Fcj)^'"^^; c,n) = 0. Furthermore, 
we can compute dx/dipc and make a change of variables with it to get 



AA[0c 



F(j/;c, n) + ]^{y~ if 



y''[2F{y-c,n)]~^'^dy (3.10) 



7 




Figure 2: A diagram of A/'[0c]- A/'[0c] and d/dcN\(t)c\ are zero at c = n. For 
c < n, there are not solitary waves. 



At present, we have only been able to evaluate (|3.10p analytically in the case 
n — 2. Using Mathematical we compute 

3c - 4 



. , , 1273 (776 + 9374 + 10572 + 35) IIT^ 

^t^^^ ^ ^^iT^T? "^V^ 

By inspection, this is strictly increasing for c > 2. The general behavior, both 
in this case and the rest, is diagrammed in Figure [51 

For n ^ 2, we justify our result with a computer plot, shown in Figure [31 
The manifold, Ai = {n^c,M[4>c\), is monotonically increasing in c for fixed n. 
This was computed by first solving F{(t)™'^^\ c,n) — Q for (f)™'^^ using Brent's 
Method. ([3T0)) was then integrated using the QUADPACK routine QAWS to 
handle the — y)^^^^ singularity. We used the GNU Scientific Library 

(GSL) implementation of these two methods, [7 . 



4 Orbital Stability and Global Existence 
4.1 Constraints of the Flow 

Under appropriate restrictions on our function space, Lc will be a positive def- 
inite operator, allowing us to conclude orbital stability. This is accomplished 
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Figure 3: A plot of log A/'[(/'c], as a function of both c and n, c > n > 1, for 
the set (c, n) e [1.5,9.5] x [1.5,9.5]. It is monotonic in both arguments, and in 
particular, increasing in c for fixed n. 
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via two constraints discussed in the following two sections. 14.1. H and 14.1.21 



4.1.1 The TV Invariant 

We would like to assume N[4)\ = A/'[0c] , and use this as a constraint, as Benjamin 
did in [3]. This is of course not true for arbitrary perturbations. Following, 
[4], we show that for cf) sufficiently close to 0c, we can find nearby solitary wave 
(j)ci such that N[(j)] = N[(t)c']- 

Lemma 4.1 Given a solitary wave <f)c, there exists (5c > such that for \\u\\^i < 
5jsf there exists a c' such that M[u -\- (pc] — M[4>c']- Furthermore, there exists a 
constant, Kc ~ Kc{4ic) such that |c — c'| < Kc . 

Proof: This is proved using the implicit function theorem. Let be the 
functional T : H\M.) x R R 

T[U, Ac] = U[cbc +U]- M[(j)c+Ac] (4.1) 

JF[0, 0] = and the Frechet derivatives are 

D.F[O,0]= (^(c-i(C-l),.), -^AA[0c]) (4.2) 

Both are bounded operators, and by Proposition [T] ,9cA/'[(/'c] 7^ 0. We may 
therefore apply the implicit function theorem for Banach spaces (see |14| for 
example), to conclude existence of an > and a mapping Q : M, 
such that if ||u||^i < 5j\r, then T[u, Q[u]] = 0. We then set c' — c + Q[u]. Since 
this map is C^, |c — c'| < K(. \\u\\jji . 

■ 

If is our perturbed (pc and ||(/)o ~ 0c|| < ^c, then we may apply Lemma 
14.11 to find pc' such that M[4)q\ = N[4)c']- Since d is a pseudo- metric on H^{R), 

A{(j)c,4'c') is time independent and 

d(0c,0c')< Uc-Pc'Wh^ <C\c-c'\<CU^-(j)c\\H^ 

which may be made arbitrarily small. So it suffices to study the stability of (p^' ■ 
Examining how M[(t)\ = M[4>c] constrains the flow, first decompose (j) as 

(f){x + Xi^,t) = (l)c{x) + u{x,t) (4.3) 

the solitary wave and a perturbation, for some xq — xo{t). Expanding J\f[(l>] = 
M[(j)c] about (j)c, 

/oo 
[/ - i^d,-)] (0c - 1) udx = ([/ - ca, (0^a,-)] i^c -i),u)^ o{\\u\\l^) 
'OO 

which we make formal in 
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any ||u||^i ^ 2' 



Lemma 4.2 Given a solitary wave (j)c, there exists a constant C such that for 
|([/-C9.(^?a,.)]('Ac-l),^i)| <C||u||^, (4.4) 



Remark 1 Since the right-hand side of (|4.4|) is quadratic in we view it 

as a near-orthogonality constraint. 

4.1.2 The Sliding Metric and the Choice of xo{t) 

As discussed in [4 with regard to the shding metric, (|2.3p . it is not true in 
general that the value of y by which one function is translated to minimize the 
norm is be finite. We will show, under some additional appropriate assumptions 
that this is the case, see [H [5] . 

Lemma 4.3 Given a solitary wave <f>c, assume that a solution to (II. ip cj), (p—l £ 
Ci([0,r) : H^{R)), satisfies 

d(0(i),0e) < (4.5) 

for t < T . Then the infimum of the function p, 

p{y) = m-)-<l^c{- + y)\\m (4.6) 

is achieved at a finite value of of y £ M. 

Proof: p{y) is obviously continuous, and because (/)(• + y) — 1 in as 
y ±00, 

lim p{y) = U - iWm + Uc - >Uc- 

y— *±oo 

But by assumption, 

inf p(y) = inf ||0(.) - + y)\\]j, = d(0,0e)' <Uc- 

So there must be some finite yo such that p(ya) < p(±oo). By the continuity of 
p, there then exists an xo{t) at each t < T such that p{—xo{t)) < p{y) for all 



Remark 2 There therefore exists a function Xq (t) such that 

d{(l),<j>,) ^ m;t) ~ Mm\m (4-7) 

Since (f)c is in fact smooth, so is p{y), and hence p'{—xo{t)) ~ 0. 

[I - dl] d.^cudx = ( [/ - dl] d,^c: u) = (4.8) 



where we used the decomposition (j4.3p . (|4.8p is a second constraint on the 
perturbation to (j)c, which together with the near- orthogonality condition (|4.4p . 
will he shown to yield local convexity of £ near (pe- 



ll 



4.2 Properties of the Linear Operator, Lc 

Here we summarize properties of and exhibit the non-positivity of L^. 

Proposition 2 (Properties of the Linear Operator Lc) The linear second order 
operator defined by (j3.5p has the following features: 

(i) Lc is self adjoint, i.e. 

Lc - Lt (4.9) 

(ii) dx(t>c is an eigenvector of Lc with eigenvalue zero 

Lcd,4>c = (4.10) 

(iii) 

Lcdccjyc = -[I~ €9. {€d.-)] {4>c - 1) - -c-i (C - 1) (4.11) 

(iv) 

{Lcdc<Pc.dc<Pc)=-^M[M (4.12) 
dc 

(v) There exists Aq < and ipo G i^(K) such that 

Lci^o = Ao^^V'o (4.13) 

i.e. ipo is the ground state of Lc with eigenvalue Ao of the generalized 
eigenvalue problem (|4.13|) . 

Proof: (i-iv) are trivial algebra and integration by parts. For (v), note that 
although Lc is not in Sturm-Liouville form, if, given h, we let h = (fi^h, 

cj^-^Lch = -cdlh + [c - ,^^"-1 - nc0r '5:. {€dx<j)c)\ h = Lch 

Lc is in standard Sturm-Liouville form, and it has a zero eigenvector, (I)cdx4'c- 
Since this has one zero crossing, by oscillation theory (see Jy^, amongst others), 
we know there exists a ground state for Lc, which we will denote by ipo with a 
negative eigenvalue, Aq. In turn, Lc has a generalized eigenvector tpo = 4'7""'Po 
with eigenvalue Aq, so 

Lcipo = Ao(/)c"V'o 

We know this is the ground state of Lc because if there existed some other tp 
with eigenvalue A < Aq, "0 = 4'"''P would be an eigenvector of Lc with eigenvalue 
A, which contradicts ipo being the ground state of Lc- 



We will now prove that with the constraints introduced in the previous 
section, p.5|) admits the estimate 

\A£\ >^{LcU,u)>C\\u\\l,+0(\\u\\%, 



12 



Defining the two vectors in L^, 

6 = [i-dl]d,(bc 



(4.14) 
(4.15) 



Proposition 3 Let 

A={feH': ||0?/IL2=1 and f±^,} 

Then 

inf (LJ, /) = 

Proof: Following. [24l [25] , let a = inf (£c/, /), taken over A. Assume a < 0, 
and let us treat this as a constrained minimization problem. From the theory 
of Lagrange multipliers, there exist J^,, such that 



(4.16) 



If Pi, is zero, then a — \q which implies is some multiple of the ground state 
V'Oi as defined in Proposition [2l But since -00 and ^ are both even functions, 

and this contradicts the assumption that /* is orthogonal to Therefore, 
/3* 0. If a = Ao, then, taking the inner product of both sides of (|4.16p with 
the ground state, 

= (A, (i, - Ao02«) = ((i, - A, V;o) = P. (6, ^o) = (0L' - 1) , ^o^ 

So Ao < a < 0. 

Let g{\) be defined as 

3(A) = ((ic-A02")'iei,Ci) 
on the interval (Aq , 0] . Note 

ff'(A) = 



^ (L,-A02")-'a,(ic-A02")-\ 



> 



so gr is increasing on this interval. Additionally, 

5(O) = -^AA[0e]<O 
by Proposition [TJ This implies that g{a) < 0. But 

5(«) = 7^(A,a) -0 



Therefore a = 0. 
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Remark 3 It is here that we see the importance of the slope condition on Af 
with respect to c. Also, we see from Figure\3l (a) and (b), and the exact com- 
putation when n — 2, that as c ^ n, dcAf[(t>c\ — > 0. There is a bifurcation point 
at c ~ n, for there are not solitary waves for c < n, but under linearization 
about (j) = 1, there are plane waves with group velocity < n; this is cartooned in 
Figure\^ (a). The sign of the derivative of this functional with respect to c was 
previously used in jl] to conclude linear instability of one dimensional solitary 
waves in two spatial dimensions of (|l.ip . when n — 3 and m = 0. 

Proposition 4 Let 

B = {feA: fU2} = {feH': 110^/11^^=1, /±e2} 
Then 

inf (L J, /) > 



Proof: Let 

a = inf (Lcf, f) 

Since ;B C a > by the previous Proposition. Assume that a = and this 
minimum is achieved at fi,. Again, by the theory of Lagrange muhiphers, 

Lcf. = + 7*6 

Taking the inner product of both ic/* with dx(t>c, 

- (A, LAfc) = {Lcf*, d.fc) = 13. {^ud^^c) + 7* (6, d^c) 

Since is an even function and dx4ic is odd, {^i,dx4>c) = 0, hence j^, — 0, 
implying 

= Ki,dx(l>c - P*dc(l>c 

Taking the inner product of f^, with , 

o = (A,a) = -/3*x-^[0c]y^o 

ac 

Therefore a > 0. 



Corollary 2 There exists a constant C > such that for all f G iJ^(M) or- 
thogonal to both ^1 and ^2 

{LcfJ)>C\\f\\l, 

which further implies 

{LcfJ)>C'\\f\\l, 
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Proof : The first part is obvious by Proposition 31 To prove the second in- 
equahty, let us express using (|3.5p as 



Then 



L, = -cd.,{^l'^d,-)+V{x) 



(LJJ) = c{cr,d,J,r,d^f) + {Vf,f) 

> mwi. \\v\\^^ ii/iii. 

> \\d,f\\l.-C-'\\V\\^^{LJJ) 



hence 



{Ljj)>{i + c-'\\v\\^^) Wd^fWl, 

and the inequahty foUows. 



Proposition 5 Given a solitary wave <j)c with c > n > 1 for (jl.ip . there exist 
positive constants C2, C3 such that for all u G satisfying 



m < t; (4.17) 



1 
2 

M[(j>c + u] = AA[0J (4.18) 

= (4.19) 



^2 given by (|4.15p . we have 



{L,u,u) > C2\\u\\l^ - CsWuWl^ (4.20) 



Proof; Let u = uj^ + M|| where 



and 



Then 



U±=U - = It - (w, ^1) ^1 



{LcU,u) = (Lc (m± + u\\) , M± + W||) = {LcUj_, u±) + 2 (icW±, + ("Ih "11) 
By Corollary [2] and Lemma [42l 



{LcU±,u±) > C = C ||w||^i - 2(m, - 2(9:jU,5a;W||) + 

\ 2 II ^ ||2 



C 



2 7-.ll_l|3 



(4.21) 
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The other two terms fohow more easily 

> -L'IImII^i (4.22) 



> -I^llull^i (4.23) 
then foUows immediately from (liT^ . I^?^ . and (li?^ . 

■ 

Remark 4 Using the estimate (I4.20p together with (|3.6p . 

> ^^2^11^1 - Qc2 + i?3) ^P- (11^*11^0 (4.24) 

T/ie polynomial p_ wiZ/ 6e important in the next section. 



4.3 Proof of Orbital Stability and Global Existence in 
Time 

Unlike non-degenerate equations, such as KdV and BBM, we need to control 
II 1/011 ^oo in our existence proof, hence the additional condition appearing in 
Theorem 12.11 It is a lack of a priori bounds on this quantity that currently 
prevents a global existence proof for general (n, m); this matter is discussed in 
[20] . However, we are able to prove, in tandem with nonlinear stability, global 
existence in time for data in a neighborhood of a solitary wave. 



4.3.1 Proof of Main Theorem 

Let a particular solitary wave (j}c be given for some c > n > 1, and let (/) be a 
solution to (|1.1[) . First we will consider the case where Af[(j)] — A/'[0c]; this will 
then be relaxed. 

For u satisfying (I4.17p . (|4.18|) . and ()4.19p . the perturbation wiU satisfy the 
two inequalities 

P-i\\u\\H,)<\A£\ <p+{\\u\\h,) 

p± defined in p.7p and (|4.24p . Since A£ is time independent, if the perturbation 
at time t = is sufficiently small, |A£| may be made arbitrarily small. Provided 
conditions (I4.17p . ()4.18p . and (|4.19p continue to hold, this will constrain ||u||^i 
through p_ . 
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Let be defined as 

which depends only on (j)c- The significance of the three quantities is: 

• 1/4 will ensure is bounded away from zero, as needed by our existence 
theory. 

• 1/2 |j(/)c — l||^i will ensure the value at which the sliding metric is mini- 
mized is finite; see Lemma 14.31 

• C2/(3C3 + 6D3) will ensure that the perturbation reamins to the left of 
the peak of the polynomial p- . 

Let e < and let (5 > be sufficiently small such that 

S <e (4.26) 
p+{6)<p-ie) (4.27) 

Letting (j){t = 0) = 0o, assume there is xo G M such that 

1100 - 0c(- + 2^0)11^1 = II^oIIhi < s 

- 1 e ([0, Tmax); H^) ■ With these choices of S, e, and e^, we will show that 
forte [0,r„,ax), d(0(t),0,) < e. 
Let 

J={t: d{(j){t'), 0c) < efor t' < t.} 
We will use X to prove the theorem by contradiction as follows: 

• Use continuous dependence upon the data of solutions of (jl.ip to prove 
that 2 is not empty. 

• Seek the maximal time Tq in X. If it is not Tmax, we will show that there 
is some time interval beyond Tq for which d{(j){t), (j)c) < 2e. 

• Prove that for any t such that d{<j){t),<j)c) < 2e, in fact d{<j){t),<j)c) < £, 
producing the contradiction. 

First we prove that X is not empty. This is an application of Theorem 12.11 

1 

00 



3 



4 

< - 
- 3 



Taking R = 3/2 ||0c — ~ "^Z^' know from part (a) of Theorem 12.11 

that ||0(i) - 1|| < 2R and ||l/0(t)||ioo < 1/a up till at least TexistX-R, ") > 

a-nd Tn^ax ^ ^cxist.- 
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By part (c) of the same Theorem, we have a constant i^cont. — ^cont.(^j 
such that 

Taking t sufficiently small, there is some time interval over which (p is within e 

of (/)c- 

Let 

To = sup J (4.28) 
Suppose To < Tmax- For any t < To, 

d(0(O,0c)<<e^ll0c-l||ffi 

Lemma [4.31 asserts there exists xo{t) £ R such that 

d(0(i), 0e) - m-,t) -M- + < £ 

Furthermore, for all t < Tq, we have again 

< ^||<^c-1|Ihi and 

We use these bounds to control how far (f> can deviate from (jjc beyond Tq. Taking 
R — 3/2 II0C — and a — 3/8, there exists Texist.(-R, a) > such that ^(t), 
for t < To, may be continued in time by the amount Texist. and will satisfy 
\\4'{t) ^ 1|| < 2_R, and ||l/0(t)||ioc < l/a. Because the solution is unique, these 
bounds apply to cf) for t < Tq + Texist.- 

With this control on norms, we let ti < Tq be a new starting point, and 
apply part (c) of Theorem 12.11 

- M- - ct + xoitMm < Ui^ti) - M- + Mti)\\m e^-' (*-*i) < ef— 
Making ti sufficiently close to Tq, we can find T2 > Tq for which 

||0(-,t)-0c(--ci + .To(ii))||^i = I|wWIIhi <2£ 

for t < T2. We claim this implies the stricter estimate <) — (pd- — ct + xo{ti))\\jji < 
e for To < t < T2. 

2e < 2ei, < ||(/>c — l||^i ; applying Lemma 14.31 to find xo{t) again, and then 
decompose our solution via (|4.3p . As noted in the remark following Lemma 
this perturbation satisfies {£,2,u) = 0. A/'[0] = JVlcpc] and 

\\u\\^, <2e<2e}- 
Therefore u satisfies the neccessary conditions to apply Proposition [5l 



F/1 



4 

< - 
- 3 
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Because < 2<j2e^, < 2C2/(3C3 + 6D3), it sits remains to the left of the 

peak of p_ ; 

\\uit)\\^, =d(0(i),0,) <e 

This holds for Tq < t < T2, Tq < T so we have a contradiction and conclude 
d((^(t),(/)c) < e for aU t < T. 

Now we relax A/'[(/)o] = A/'[0c]- First apply Lemma WJ\ to </)c, to find 5c that 
will define the neighborhood about the origin where the initial perturbation 
must reside. Let K be the constant such that 

Uc-M\h^ <K\c-c;\<kKc\\u\\H,=K\\u\\H, , 

where Kc is from Lemma |4. II We will seek a new wave (f)c' with which to apply 
the preceding argument; however, Si, and the polynomials p± will be determined 
by 4>c' ■ But (f)c' is determined by 0o and we have not yet found all the bounds 
(jjQ must satisfy. Uniform control is needed. Using Q, the implicit function 
associated with Lemma HTTl let 

Cmin = min g[u] 

Cmax = max g[u] 

\\u\\^l<Sa 

and let 

el = min £*((?!'c) 

CG[Cmin,Cniax] 

Similarly, the coefRcients of p± may be chosen such that 

P-(hllHi) < |f[</.c + M]-f[0c]| = |Af| <pV(||m||^0 

for all c £ [cmin , Cmax] , 11^11^1 < min{(5c, ^} (the i is required by Lemma [??T|) . 

Given (j>c, let e'^ be as above. Let e^, = el(l + K). Take e < e*, and then set 
e' = e/(l + -RT). Choose 5 such that 

d < Sc (4.29) 
(l + i^)(5 < e' (4.30) 
pV((1+X)5) < p'_{e') (4.31) 

and let 5' = (5(1 + i^). 
Assume 

1100 - 0c(- + xo)\\fji < 6 

Let 0c' be the neighboring solitary wave for which A/'[0o] — ■^['Pc']- Note (po is 
also close to 0c' , 

||0o(-)-0c'(- + a;o)|lHi < ||0o(-)-0c(- + a;o)||^i + ||0c' -0c||ffi < (l+X)-? = 5' 

Because <5' < e' and p+((5') < p'_{e'), e' < e'^, we may apply the previous 
argument, to conclude 

d(0(i),0c') <e' 
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forte [0,r„,ax). Finally, 
dm),(t>c) < d(0(i), 0c') + d(0e', -^c) <'e + KS<',+ KS' < (1 + K)e' 



4.3.2 Proof of Corollary E 

Now we will prove global existence in time for a data in a neighborhood of a 
solitary wave. Given a solitary wave 0c, e < £*, and any function such that 

Il0o(-) - <Pc{- - Xo)\\fji < S 

Let Tmax be the maximal time of existence of 0, the solution emanating from 

00. Suppose Tinax < oo. 

By Theorem for all t < T^ax, 



d(0(t),0c) <e 



which implies 



4){t) > 1 - e > 1/2 
- 111^1 < e + - lllffi < 2 - lll^i 



and hence 



< 



+ 2 



for all t < Tmax. But 6 is in time, and since Tmax < 



lim 



a contradiction. Therefore Tmax = oo, and furthermore, by Theorem 12.41 again. 

d(0(t),0,) <e 

for all time. 



5 Relation to Compacton Equations 

Compactons, robust compactly supported solitary waves, were first identified in 
[IB] in a generalization of KdV, K{m,n), 



dtu + (w™) + ^dl K) = 
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y 

Figure 4: Compacton solutions to with n — 2, m = —2 for c = 5 and 

c = 8. Faster compactons are both broader and taller. 



K^m, n) has been further generalized to Ci(to, o + 6), 

dtu + 9, (m™) + [m'^c'^ = (5.1) 

discussed in [TS] , which are known admit traveling wave solutions of speed c > 
0, compactly supported in space. In turn, this has been generalized to the 
multidimensional case, Cjv(to, a + b), in p!5l [T7] . 

Like (jl.ip . (|5.ip is a nonlinearly dispersive wave equation. Indeed, (|l.ip may 
be written as 

+ 9. m - T^d, [^""d^dt (<^i-")] = (5.2) 
1 — m 

which, upon setting a — n and 6 = 1 — m, resembles ()5.ip up to one dx becoming 
—dt- It was the recognition of this similarity between these two equations that 
led the authors to discover that just as (|5.ip is Hamiltonian for 5 — a = 1, so 
too is (|l.ip when I — m — n — 1, n + m — 0. (jl.ip may be interpreted as a 
generalization of the BBM equation, just as (|5.ip is a generalization of KdV. In 
addition, (|l.ip also possesses compacton solutions, pictured in Figured] 

Equation (|5.ip has a second set of Hamiltonian cases when b — a = ^ , which 
is related to the property that the density of the generalized momentum of (|5.ip 
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is a mapping of a (j5.ip into another such equation. The generaHzed momentum 
of (|3.ip . not being a monomial, lacks this property. Note also that 

could be further generalized by letting the the nonlinearity in the (<^")x vary 
independently of the dispersive term; however, for physical reasons related to 
its derivation, we leave it as is. 

We now construct compactons. Starting with (|5.2p . we integrate it up, as- 
suming there exists j/max > 0, such that for \y\ > i/max, Qc = 0. Then 

- cQc + Q" + = for |y| < 2/„,ax (5.3) 

1 — 771 " ^ ' 

with boundary conditions Qc ~ a.t \y\ — ^max- 
Introducing the scaling 

Q^{y) = c^[/(Vl - mc^^^y)^ = c-^U{^)^ (5.4) 
into (|5.3p . [/ solves 

- [/-^ + 1 + a|C/ = (5.5) 

Remark 5 Assuming n > 1, the scaling, (|5.4p . highlights three distinct regimes 

• If n > m, faster waves are taller and narrower. 

• If n < m, faster waves are taller and broader. 

• If n — m, faster waves are taller, but all waves have the same width. 



Integrating (|5.5p again, 



where the potential, V, 



^ ''^^ \ U-\ogU, if/3=l 

and /3 = (n — 1)/(1 — rri). 

Representative cases for/3>l,0</3<l and f3 < are illustrated Figure 
[51 Compactons are homoclinic orbits in the {U, U') phase plane connecting the 
equilibrium point (0, 0) to itself in finite time. This corresponds to V having a 
potential well between U — and some C/max > 0. Such a well exists only for 
< /? < 1; see Figure[5] 

In terms of n and tti, for 7t,>1, 0</3<l corresponds to ti + 777, < 2, which 
obviously includes the Hamiltonian case of l|1.2p . (|1.3p . (II. 4p : 77 + 777, = 0. With 
regard to Remark [51 because we have assumed 77 > 1, we are in regime where 
faster waves are both taller and broader. 
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/3>1 




^ ^ ^ 1 




/3 < 


\ 




, . • ■ ^ 







Figure 5: The different regimes of P for different /?. We exclude = because 
P(/3 = 0) = 0. 



We note that several examples of compactly supported solutions of (jl.ip 
were previously examined in [23, 22J. These included the exponents (n,m) = 
(3,0), (4,0). The authors argued against the realization in nature of such solu- 
tions due to a stress singularity. 

Some solutions of (|5.5|) . given implicitly by 



24J7 + 10J74/5 + 50C/3/5 - 125[/2/5 ^ ^ ^ , ^. 
= t for n — i(5.6) 

= lel for P = i(5.7) 



V3[/2/3 _ 2U 

^/U {2^J^/U~2\og (^VVU - 2 + + U^^^ - 2VU 
^^-U/2 

and 

C/e (o. 



lel for/3=i(5.8) 



1-/3 



The compactons are a type of weak solution to but the exact notion is 

still imprecise. Near the left edge of the compacton, [/ ^ |^ + ^maxl^^^^'*''^^ H{^+ 
■Cmax), H{x) the heavyside function. In the Hamiltonian case Qc ~ \y + Vmaxl^^'^ H{y+ 
J/max)- For n > 2, this will not have a square integrable derivative. They do 
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satisfy the following definition, previously given in for solutions of (|l.ip that 
go to zero. 

Definition 5.1 (f){x,t) is a solution of p.ip if 



oo /'OO 



-dt^ix, t)(t){x, t) - dxi'ix, t)(t>{x, ty 

^ dxil^{x,t)(t){x,t)''dxdt{(t>{x,tf-''Mdxdt = Q (5.9) 



1 — TO 



for all ipix,t) g C^(R x K+), and such that both 

(f>{x,t) and (j){x,t)''dxdt{<p{x,ty-"') 
are in Lj^^{M.) in the x coordinate. 

Proposition 6 Qc is a solution of (jl.ip in the sense of Definition \5.1 



Proof: For \x — ct\ > j/max, Qc = 0, and satisfies (|5.3p pointwise. For \x — ct\ < 
2/inax, Qc is smooth and solves (|5.3p in the classical sense. Now consider y = x—ct 
near -?;„iax- In this neighborhood, Q^a^ (Qj-™) ~ |y + ymaxl^^^""™'' ^^(y + 
2/max)- So for 71 — TO > 0, this is a continuous function and (|5.3p also holds 
pointwise. Therefore, 



oo /"OO 



-dtip{x,t)Qc{x - ct) ~ dxi^{x,t)Qc{x - ct)"^ 

t)Q,{x - ct^dl (Q,(x - ci)!-") 



1 — TO 

oc i*oo 



dxdt 



{—dtip{x, t) — cdxip{x, t)) Qc{x — ct)dxdt 



But this is the weak form of the transport equation dfU + cdxU = 0, which Q^. 
solves. So the integral is zero for all test functions and Qc is a weak solution in 
this sense. 



Finally, if instead of having the Hamiltonian as in (|1.3p , we set 

1 



H 



n+l 



dx 



and replace the generalized momentum M with 



N ■ 



0'" [dx^Y + dx 



then these formally conserved quantities are finite for the compactons. The 
compactons are then critical points of the energy functional £ = TL + cj\f. This 
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suggests the possibility of an analogous stability argument as for the solitary 
waves. One can take the second variation and formulate the spectral problem 



just as in Section [5T^ However, a well-posedness theorem for solutions of (|l.ip 
that vanish outside a compact set must be formulated before this is pursued. 

6 Remarks and Open Questions 

We have presented a new class of Hamiltonian PDEs with orbitally stable soli- 
tary waves. A consequence of our stability analysis is the extension of well- 
posedness results for this system in ^D] to a neighborhood of any solitary wave. 
As noted, except for the case n = 2, this result is currently only valid up to 
the acceptance of numerical computation and estimation of the slope of the in- 
variant A/'[(/)c]- We also observed that our equations have compacton solutions. 
These are solitary traveling waves, whose spatial support is compact. We show 
that these compactons solve the evolution equation in a weak sense. 

Compactons warrant further examination. Formally, compactons are critical 
points of the functional £[(j>], defined in section [Sj A well-posedness theory in a 
function space, with respect to which the mapping (j) i-^ £[(j)] is continuous, and 
a spectral analysis of the second variation about a compacton analogous to that 
for the solitary waves of section [U which would imply stability of compactons. 
This is an interesting open problem. 
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A Generalized Momentum 

The Lagrangian density, C, for the Hamiltonian case of (|l.ip is given by ^px = 



2 n -|- 1 2 



(jl.ip is then Euler-Lagrange equation corresponding to S J J Cdxdt. 




(A.2) 
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and 

the generalized momentum, p.ip . of (jl.ip in the case n + m = 0. 
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